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1 . 0  SUMMARY 


The  overall  preliminary  theoretical  flutter  analysis  of  the  Ryan 
XV- 5A  Lift  Fan  Aircraft  is  presented  in  three  volumes.  This  is 
Volume  II,  which  reports  the  investigation  of  the  flutter  charac¬ 
teristics  of  the  empennage.  Volumes  I  and  III  detail  the  analysis 
of  the  wing  and  control  surfaces,  respectively. 

Results  of  the  empennage  flutter  analysis  indicate  a  possible 
flutter  situation  existing  within  the  design  flight  envelope  of 
the  XV- 5A  aircraft.  The  flutter  mode  is  characterized  by 
horizontal  tail  pitch  coupled  with  basically  symmetric  horizontal 
tail  bending  for  low  values  of  pitch  frequencies  (40  cps)  and  for 
pitch  coupled  with  horizontal  tail  torsion  at  pitch  frequencies 
greater  than  40  cps.  The  required  pitch  frequency  would  be  of 
the  order  of  approximately  55  cps  for  a  flutter  velocity  margin 
of  15%  of  the  design  limit  dive  speed  at  sea-level,  which  is 
the  critical  altitude.  Establishment  of  the  actual  pitch  frequency 
through  experimental  methods  is  required  for  definite  determination 
of  flutter  margins.  The  analyses  presented  in  this  report  should 
be  taken  as  a  preliminary  evaluation  of  the  empennage  of  the  XV- 5A 
aircraft.  Subsequent  analytical  and  experimental  investigations 
should  then  be  directed  toward  finalization  of  the  actual  flutter 
characteristics  of  the  XV-5A  aircraft,  with  due  regard  for  the 
findings  of  the  studies  presented  in  this  report. 


1 


'IT 


2.0  INTRODUCTION 


This  report  presents  the  results  of  the  preliminary  empennage  flutter 
analysis  of  the  U.  8.  Amy  XV-5A  lift  Fan  Researoh  Aircraft,  illustrated 
in  Figure  1. 0  of  Reference  1.  The  XV- 5 A  is  a  V/8TOL  airoraft  designed 
for  researoh  flight  testing  of  the  General  Electric  X-353-8  Lift  Fan 
Propulsion  System. 


The  XV- 5 A  aircraft  is  characterized  by  the  existence  of  a  T-tail  with 
manual  type,  conventional  aerodynamic  control  surfaoes  ( elevator  and 
rudder).  Suoh  tails  are  prone  to  exhibit  flutter  instabilities,  unless 
extreme  oare  is  taken  during  the  design  stages  to  insure  a  flutter-free 
structure.  Further  difficulties  are  presented  due  to  the  fact  that  the 
horizontal  tail  is  pivoted  to  provide  for  longitudinal  trim  and  is  powered 
by  an  irreversible  Jackscrew  with  power  applied  to  the  Jaoksorew  through 
hydraulic  oircults. 

The  analysis  of  the  empennage  was  facilitated  by  the  use  of  simple 
beam  Idealization  of  the  horizontal  and  vertical  stabilizers  and  fuse¬ 
lage.  Aerodynamic  representation  was  accomplished  through  the 
use  of  two-dimensional  strip  theory.  Parametrto  investigations  of 
important  items  suoh  as  the  horizontal  tail  pitch  frequency,  elevator  and 
rudder  rotational  frequency,  and  horizontal  tail  rigid  body  yaw  and  roll 
were  carried  out  for  establishment  of  oritioal  areas. 


3.0  METHOD  OF  APPROACH 


The  method  of  approach  followed  in  the  flutter  analysis  of  the  empennage 
is  that  of  the  Rayleigh-Ritz  approach,  or  what  might  be  termed  the  modal 
approach  to  the  flutter  problem.  In  this  approach  the  equations  of  motion 
are  formulated  in  terms  of  carefully  chosen  modal  functions.  The  modal 
functions  chosen  to  represent  the  deformation  of  the  structure  are  the 
calculated  normal  modes  of  vibration  of  the  particular  component,  i.e. , 
horizontal  and  vertical  stabilizers  and  fuselage. 


3. 1  LAGRANGE'S  EQUATIONS  OF  MOTION 


In  applying  the  modal  approach  to  the  flutter  problem,  Lagrange's  equa¬ 
tions  of  motion  are  convenient  for  obtaining  the  equations  of  motion 
necessary  for  a  flutter  solution.  In  general,  Lagrange's  equations  may 
be  expressed  as, 


d  d  T  |  9U  8D 

dt  dq.  9q.  8q.  j 

J  J  J 


where: 


th 

4^  is  the  j  generalized  coordinate 
T  is  the  kinetic  energy  of  the  system 


(1) 


U  is  the  potential  energy  (internal  strain  energy)  of  the  system 

D  is  a  dissipation  function  and  is  composed  of  damping  forces  in  the 
structure 

Q  is  the  generalized  force  applied  to  the  system  in  the  coordinate. 

The  kinetic  energy  expression  T  may  be  written  as, 

T  =  Y  q  '  M  q  (2) 

where  M  is  termed  the  generalized  mass  matrix  and  is  composed  of 
terms  which  are  products  of  the  lumped  masses  and/or  mass-inertias 
and  the  chosen  modal  functions. 


The  potential  energy  expression  U  may  be  written  as 

u  -  |  q'Ea 

where  K  is  termed  the  generalized  stiffness  matrix  and  is  composed  of 
terms  which  are  products  of  the  frequency  of  oscillation  squared  of  a 
particular  mode  and  its  appropriate  generalized  mass. 

The  dissipation  function  D  is  due  to  the  damping  forces  acting  in  the 
structure.  These  may  be  considered  as  forces  in  phase  with  the  velocity 
and  proportional  to  the  restoring  forces  acting  in  the  structure.  There¬ 
fore, 


« 

(3) 


D  =  - 


_g 

2u) 


S'Kq 


(4) 


where  g  is  the  coefficient  of  structural  damping  and  w  is  the  frequency  of 
oscillation. 


The  generalized  force  Q,  is  due  to  the  oscillatory  aerodynamics  acting  on 
the  structure  and,  for  all  force  components  j,  Qj  may  be  written  as, 

q  =  co 2  A  q  (5) 

where  A  is  called  the  generalized  aerodynamic  matrix. 

Applying  Equation  l  to  Equations  2,  3,  4,  and  5,  and  assuming  harmonic 
moation  (q  =  qGeiwt),  the  equations  of  motion  for  a  flutter  solution  are, 

2  —  —  2 
-0)  M  q  +  K  q  +  ig  K  q  =  0)  A  q 

or 

-  -%  Ea  =  -  A.a 

or 

(i  +ig)  S )  a  =  o 


Ma  -  -4  £a 


Rearranging, 


M  +  A - - 

'  “  CD2 


6 


i 


or 

(M+A-ftK)q=o  (6) 

/  2 

where  ft  =  l/u>  (1  +  ig). 

For  a  solution  to  exist,  the  determinant  of  the  matrix  of  Equation  6  must 
be  equal  to  zero,  i.e., 

|M  +  A-  ftK|=o  (7) 

Solution  of  Equation  7  by  standard  means  yields  the  required  complex 
eigenvalues  from  which  the  flutter  speed,  frequency  and  damping  values 
are  determined. 

3.1.1  Generalized  Mass  Matrix 

This  section  is  devoted  to  the  derivation  of  the  generalized  mass  matrix 
in  a  general  sense  to  which  particulars  are  applied  in  subsequent  sections. 

Assume  a  given  element  of  mass  dm,  as  shown  in  Figure  1,  which  repre¬ 
sents  a  portion  of  the  i^  section  of  a  lifting  surface.  A  rectangular  car¬ 
tesian  coordinate  system  is  fixed  to  the  surface  at  point  i,  which 
represents  a  point  on  the  elastic  axis  of  the  surface. 


__  .  Section  i 

Point  i  > i 


Z,  w 

Figure  1  i^  Section  -  Lifting  Surface 
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The  displacement  of  the  element  dm,  in  terms  of  displacements  of  point  i, 
may  be  written  as, 

u  =  <$x  -  6^z  -  6zy  (8) 


v  =  <5  -6  z  +  (5  x 

Y  X  Z 


*  =  v  V+  V 


Equations  8,  9,  and  10  may  be  written  in  matrix  form  as 


1.0  0  0  0  -z  -y 


0  1.0  0  -z  0  x 


0  0  1.0  y  x  0 


A  =  B  6' 


The  kinetic  energy  of  the  i  section  is 


a  •  •  I  • 


2T  = 


L  u  v  wj  u  I  dm 


Section 


And  by  substituting  Equation  11, 


2T  = 


-/ 


1(1)'  „  1(1)  . 
d  B1  B  d  dm 


Section 


=  i(1)'  / 


B'  B  dm  d 


Section 


8 


Lift*  *c 


Multiplying  out  the  matrix  product  B'  B  yields, 


1.0 

0 

0 

0 

-z 

-y 

0 

1.0 

0 

-z 

0 

X 

0 

0 

1.0 

y 

X 

0 

0 

-z 

y 

2  x  2 

y  +  Z 

xy 

-xz 

-z 

0 

X 

xy 

2  2 
x  +  z 

yz 

-y 

X 

0 

-xz 

yz 

2 

x  +  y 

Therefore,  the  kinetic  energy  expression  may  be  written  as 


2T  =  A(i)y 


Section 


1.0  0  0  0 

0  1.0  0  -z 


-z 


-y 


0 

0 

1.0 

y 

X 

0 

0 

-z 

y 

2  2 
y  +z 

xy 

-xz 

-z 

0 

X 

xy 

2  2 

X  +z 

yz 

-y 

X 

0 

-xz 

yz 

2  : 
x  +y 

dm  «<*> 


Integrating  over  the  itXl  section  results  in  the  following  expression  for  the 
kinetic  energy  of  the  lumped  parameter  system  of  the  ixXl  section, 


2T 


6<1)’ 


m 

0 

0 


0 


m 


-Syz  -SZy 


-s 


xz 


0 


sz. 


m  S 


XY 


0  -SXZ  SXy  XX 


-s 


YZ 


-Sry  Sr, 

Zy  ZX 


s 


YX  PXY 

0  _pxz 


SYX  0 
^XY  _PXZ 

PYZ  !Z 
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or 


2T  =  «(l>'  J(1)  6(1) 

where 

m  is  the  mass  of  section  i  Lumped  at  station  i 


(14) 


Z 


is  the  mass  moment  of  inertia  of  section  i  about 
the  X,  Y,  Z  axes,  respectively 


P^Y»  ^xz*  ^YZ  *8  ma8S  Pro<Juct  °f  inertia  of  section  i  about 

the  X-Y,  X-Z,  and  Y-Z  axes,  respectively 


is  the  static  mass  moment  of  section  i  about  the 
X-axis,  with  distances  measured  along  the  Y  and 
Z  axes,  respectively 


is  the  static  mass  moment  of  section  i  about  the 
Y-axis,  with  distances  measured  along  the  X  and 
Z  axes,  respectively 


S 


is  the  static  mass  moment  of  section  i  about  the 
Z-axis,  with  distances  measured  along  the  X  and 
Y  axes,  respectively. 


Normally,  the  surface  to  be  analyzed  is  swept;  therefore  a  coordinate 
transformation  is  required,  since  the  analysis  is  based  upon  elastic- 
axis  displacements.  Consider  the  coordinate  system  shown  in  Figure  2, 
where  ♦g  a  is  the  sweep  back  angle  of  the  elastic-axis  of  the  surface. 

* 


10 


The  displacements  6^  given  in  Equation  14  may  be  written  in  terms  of 
elastic-axis  displacements  as  follows: 


Substituting  Equation  15  into  Equation  14  and  performing  the  required 
matrix  multiplication, 

2T  =  6^'  T'  T  6^ 

or 


11 


4 


Since 


1 

u 


2 

cos 


* 

e.a. 


2 

+  I,  sin  4» 

Y  e.a. 


+  2  1>  r  sin  *  cos  >1/ 

XY  e.a.  e.a. 


I 

v 


2 

=  I  sin  * 

X  e.a. 


_  2 
+  I  ,  cos  ♦ 

Y  e.a. 


-2  P,r,rsin>|/  cos\|/ 

XY  e.a.  e.a. 


P 

uv 


(L-L,)sin4»  cos* 
Y  X  e.a.  e . 


U  • 


+  PXY<C0S\a. 


-  sin  ) 
e.a. 


S 

u 


COB* 

e.a. 


sin* 

e.a. 


S 

v 


cos  S' 
e.a. 


sin  4* 
e.a. 


Equation  16  may  be  written  as 

o 


2T  = 


m 


0  -S  sin*  -S  cos* 

Y  e.a.  Y  e.a. 

Ct  *-* 


-s 


z. 


m  0  -S  cos*  S„  sin* 

Xz  e.a.  Xz  e.a. 


m 


•8  Bin*  -S  cos*  8 
Y  e.a.  X  e.a.  u 


-8  COB*  S„  sin*  8 

Y  e.a.  X  e.a.  v 
z  z 


uv 


uv 


-PXZCOB*e.a. 

+PYZBln*e.a. 


P  sin* 

XZ  e.a. 

+P.  COB* 

YZ  e.a. 


-8. 


~P  cob*  Pv,.8ln  + 

XZ  e.a.  XZ  e.a. 

5„  0  +P„„«ln*  +P„,  cos* 

Z  YZ  e.a.  YZ  e.a. 

A 


z 


;<2) 


Since  we  are  dealing  with  surfaces  lying  mainly  in  the  X-Y  plane  (X-Z 
plane  in  the  case  of  the  vertical  stabilizer),  all  distances  in  the  Z  direc¬ 
tion  are  small  compared  to  X  and  Y  distances,  and  the  mass  properties 
involving  these  distances  may  therefore  be  set  equal  to  zero. 


% 
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i 


Equation  17  therefore  reduces  to  the  following: 


With  the  mass  distribution  taken  as  in  Figure  2,  i.e.,  streamwise  strips, 
the  following  may  be  substituted  for  the  static  moment  off-diagonal  terms 
in  Equation  18: 


S„  =  S  sin  +  S  cos  * 

Z  u  e.a.  v  e.a. 

A 


Equation  18  therefore  takes  the  following  form: 


13 


or 


2T  -  i(2>’  J<2>  6<2) 


(19) 


For  mass  distributions  taken  perpendicular  or  nearly  perpendicular  to  the 
elastic-axis,  various  terms  of  the  mass  matrix  J^2)  take  on  the  following 
values: 


S  =  0 
u 


P  =  0 
uv 


S 


-S  sin  * 
v  e.a. 


S„  =  S  cos  ’P 
Z  v  e.a. 

A 


Equation  18  therefore  reduces  to  the  following  form: 


•  (2)' 
2T  =  6'  ' 


m 


0 

0 

0 


S  sin* 
v  e.a. 


0 

m 

0 

0 

0 


S  cos\p 
v  e.a. 


0  0 


m  0 

0  I 
u 

S  0 

v 

0  0 


0 

0 

s 

v 

0 

I 

v 

0 


S  simp 
v  e.a. 

S  cosvp 
v  e.a. 

0 

0 

0 


;<2) 


or 


2T  =  F  i<3»  i<2> 


(20) 


Equations  19  and  20  are  written  for  motion  of  the  basic  surface;  for  in¬ 
clusion  of  control  surface  degrees  of  freedom,  additional  displacements 
must  be  included  and  the  corresponding  elements  added  to  the  mass 
matrix.  (Mass  and  mass-inertia  effects  of  a  control  surface  are  included 
in  the  basic  surface  terms  of  Equations  19  and  20,  since,  in  effect,  the 
control  surface  in  a  locked  position  is  part  of  the  basic  surface.) 


14 


Let  6W  be  the  angular  deflection  of  the  control  surface  about  its  hinge- 
c 

line,  relative  to  the  basic  surface  as  shown  in  Figure  3, 


STATION  i 


Figure  3  Control  Surface  Coordinate  System 

where  the  w-axis  is  the  hinge-line  of  the  control  surface  and  is  non¬ 
parallel  to  the  elastic-axis  as  shown. 

Upon  including  the  control  surface  degree  of  freedom  6  , 

takes  the  following  form: 


Equation  19 


ZT 


2T-I  6  6  6  6  6  6  6  I 

L  X  Y  Z  u  v  Z  w„  J 


in 

0 

0 

0 

0 

0 

0 

8  Bin* 

U  ti.  tt. 

0 

in 

0 

U 

0 

+8  cob* 

0 

6V 

v  e.  u. 

Y 

0 

0 

ni 

8 

8 

0 

8 

6 

u 

V 

w 

/ 

V 

0 

u 

8 

I 

p 

0 

P 

6 

u 

u 

uv 

uw 

r 

u 

0 

0 

S 

I1 

1 

0 

P 

6 

V 

uv 

V 

vw 

V 

V 

8  Bin* 

U  ti .  J  . 

0 

*8  COB* 

0 

0 

0 

0 

6,. 

v  e.a 

z 

/ 

0 

0 

8 

P 

P 

0 

1 

L 

w 

uw 

vw 

w 

w 

_ 

c 

c 

c 

e 

c 

or 


2T  =  a<2>'  j«2>  a<2> 

And  Equation  20  becomes 


(21) 


\  6  6  i  6  6  6  t  I 

[_  x  y  /  u  v  /.  wu  J 


— 

— 

Ill 

0 

0 

0 

0 

8  Bln* 

v  t*.a. 

0 

0 

III 

0 

0 

0 

8  cob  4* 
v  ti.  a. 

0 

0 

u 

m 

0 

8 

0 

8 

6 

V 

w 

Z 

e 

0 

0 

U 

I 

0 

0 

1* 

u 

uw 

c 

u 

0 

0 

8 

0 

1 

0 

p 

6 

V 

V 

vw 

V 

c 

8  *11)4* 

8  cum 4* 

0 

0 

0 

\„ 

0 

6 

a 

9 

> 

v  c.a. 

z 

Z 

0 

0 

8 

p 

p 

0 

1 

A 

w 

uw 

vw 

w 

w 

c 

c 

c 

c 

c 

or 


2T  =  6(2)’  J<3>  i(2) 


(22) 


where: 

I  is  the  mass  moment  of  inertia  of  the  control  sur- 

w 

c  lace  portion  of  section  i  about  the  w  axis  (control 

surface  hinge-line) 
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p 

uw 

c 


is  the  mass  product  of  inertia  of  the  control  surface 

portion  of  section  i  about  the  u  and  w  axes 

(  =  1  sin  <#  "r  P  cos<p->-S  i  sin  4*  ) 

w  sw  w  e.a. 

c  c  c 

is  the  mass  product  of  inertia  of  the  control  surface  por¬ 
tion  of  section  i  about  the  v  and  w  axes 

(  =  I  cos  q>  -  P  sin  <p  +  S  i  cos  4*  ) 

w  sw  w  e.a. 

c  c  c 

is  the  static  mass  moment  of  the  control  surface  portion 
of  section  i  about  the  w  axis  (control  surface  hinge-line) 


P 

sw 

c 


is  the  mass  product  of  inertia  of  the  control  surface 
portion  of  section  i  about  the  s  and  w  axes 


<0  =  'I' 

e.a. 


* 

,  where  is  the  sweep  back  angle  of  the  control  sur¬ 
face  hinge-line 


I 


is  the  streamwise  distance  from  the  v-axis  to  the  w-axis 
at  station  i 


Having  established  the  kinetic  energy  in  terms  of  elastic  axis  deflections, 
we  may  now  write  these  deflections  in  terms  of  generalized  coordinates 
as 


6^  =  <£  (23) 

where  the  <p' s  are  the  chosen  modal  functions  and  the  matrix  </>may  be 
termed  the  generalized  displacement  matrix. 

Substituting  Equation  23  into  Equations  21  and  22  and  summing  over  all 
the  sections  of  the  surface,  the  kinetic  energy  expression  for  the  surface 
is  then, 

(2).- 

2T  =  q’  (p1  JV  (p  q  =  qr  M  q  (24) 


or 


(3) 

2T  =  q*  (p'  «T  <p  q  =  q1  M  q  (25) 

where  the  matrix  M  is  the  generalized  mass  and  the  mass  matrix  takes 
the  form  as  in  Equation  21  or  22  but  as  diagonal  matrices  for  each  mass 


parameter.  The  generalized  displacemem  matrix  ^is  now  of  the  follow¬ 
ing  form, 


<P  = 


”11 

^21 

’’SI 


<p 


41 


<P 


51 


^61 


*71 


where  the  elements  are  column  matrices  relating  the  displacements 

of  the  elastic  axis  /  <5  ,  6  ,  6  ,  etc.\  to  the  "n"  generalized  coordinates. 

V  X  Y  Z  / 

3.1.2  Generalized  Stiffness  Matrix 


In  any  linear  elastic  system,  the  forces  acting  may  be  related  to  the 
deflections  by  the  following  equation; 

F  =  K  6  (26) 

where  the  matrix  K  is  the  stiffness  influence  coefficient  matrix,  an  ele¬ 
ment  of  which  yields  the  force  at  point  i  due  to  a  unit  deflection  at  point  j. 

The  total  internal  strain  energy  of  the  system  therefore  may  be  written 
as 


U  =  ±  6'  K  6  (27) 

Expressing  these  deflections  in  terms  of  generalized  coordinates,  as 

6  =  <pq  (28) 
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iMM 


k 


and  substituting  into  Equation  27,  we  obtain  an  expression  for  the  inter¬ 
nal  strain  energy  which  is  in  terms  of  the  generalized  coordinates  of  the 
system : 

U  =  -  q'  &  K  ipq 

=  |q'  Kq  (29) 


where  the  matrix  K  is  the  generalized  stiffness  matrix. 

Equation  1,  Section  3.1,  states  that  there  is  one  Lagrange's  equation  of 
free  undamped  vibration  for  each  generalized  coordinate  (j=l,2,3  —  n). 

Therefore,  for  the  coordinate, 


d  9T  8U 

—  —  +  -  =  Q 

dt  9q.  8  q.  j 

J  J 


(30) 


Applying  Equation  30  to  Equation  24,  Paragraph  3. 1. 1,  and  Equation  29, 
Paragraph  3.1.2,  results  in  the  following:  (For  the  generalized  force  Q 
equal  to  zero) 


J 


1  M.,  M.0  —  M.  1 
L  Jl  j2  jn _ 1 

f 

\ 

+  LKjl  KJ2  -  KjnJ 

r  i 

qi 

< 

'q2 

1 

1 

1 

1 

St 

J 

'  i 

q2 

1 

1 

| 

1 

1 

l 

Im 

(31) 


If  we  assume  all  coordinates  except  the  equal  to  zero,  the  system 
moves  in  one  coordinate  and  has  a  natural  frequency  appropriate  to  that 
coordinate.  This  is  usually  called  the  uncoupled  frequency  of  the 
coordinate. 


Therefore,  Equation  31  reduces  to 


M..  q.  +  K  q. 
JJ  J  jJ  J 


0 


If  we  assume  harmonic  motion  in  the  system, 
2 

q.  =  -w.  q. 

J  J  J 


(32) 


(33) 
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Substituting  Equation  33  into  Equation  32, 


2 

-GJ 

j 


M..  q. 

jj  J 


+  K..  q.  =  0 
JJ  J 


(co .  M..  -  K. . ) q 
J  JJ  JJ/  J 

Since  q^  ^  0, 


M.. 

JJ 


0 


(34) 


(35) 


And  the  generalized  stiffness  matrix  K  reduces  to  the  following  diagonal 
form: 


Equation  35  shows  that  zero  elastic  coupling  exists  between  the  coordi¬ 
nates  as  is  the  case  with  uncoupled  or  orthogonal  modes. 


3.1.3  Generalized  Aerodynamic  Matrix 


Consider  the  following  displacements  of  the  i^  strip  shown  in  Figure  4 
which  represents  a  streamwise  portion  of  the  lifting  surface  as  shown: 


6 

u 


6 

v 

6 

Z 


Deflection  of  the  elastic-axis  in  the  X-direction 

Deflection  of  the  elastic-axis  in  the  Y-direction 

Deflection  of  the  elastic-axis  in  the  Z-direction 

Angular  rotation  of  the  surface  about  the  u-axis 

Angular  rotation  of  the  surface  about  the  v-axis 

Angular  rotation  of  the  surface  about  the  Z-axis 

Angular  rotation  of  the  control  surface  about  the  w-axis, 
relative  to  the  basic  surface 
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i 


X  u 


Figure  4  Ith  Strip  Aerodynamic  Simulation 

These  displacements  may  be  expressed  in  terms  of  the  generalized  coor¬ 
dinates  as 

(2) 

6'  =  <£  a  (36) 

Now,  since  the  above  displacements  represent  displacements  in  terms  of 
elastic-axis  displacements,  we  may  transform  them  so  as  to  reflect  dis¬ 
placements  about  axes  perpendicular  and  parallel  to  the  free  stream 
through  the  i^1  strip. 

Let 

6^  =  Deflection  of  the  elastic-axis  in  the  X  direction 

6y  =  Deflection  of  the  elastic-axis  in  the  Y  direction 

6  =  Deflection  of  the  elastic-axiB  in  the  Z  direction 

Z 

6  =  Angular  rotation  of  the  surface  about  the  X-axis 

A 

dy  =  Angular  rotation  of  the  surface  about  the  Y-axis 

6  2  =  Angular  rotation  of  the  surface  about  the  Z-axis 

6y  =  Angular  rotation  of  the  control  surface  about  the  Y-axis, 
c  relative  to  the  basic  surface 


J  •**...  i 


Therefore, 


6 


(1) 


T<i>  a<2> 


(37) 


where  the  matrix  T  is  a  transformation  matrix  relating  displacements 
about  the  u,  v  axes  to  displacements  about  theX,  Y  axes  through  the  same 
point. 


Since  the  aerodynamic  forces  are  written  for  quarter-chord  displace¬ 
ments,  a  transformation  of  displacements  from  the  elastic-axis  to  the 
quarter-chord  is  required.  Accordingly,  let 

h  =  Vertical  displacement  of  the  quarter-chord  point,  plus  positive 
down 

a  =  Angular  rotation  of  the  surface  about  a  line  parallel  to  the  Y-axis 
through  the  quarter-chord  point,  plus  leading  edge  up 

ft  -  Angular  rotation  of  the  control  surface  about  a  line  parallel  to  the 
Y-axis  through  the  control  surface  hinge-line,  relative  to  the  basic 
surface,  plus  trailing  edge  down, 

And  the  transformation  is 


h 

O' 

ft 


t<3>  6<1) 


(38) 


The  aerodynamic  forces  acting  per  unit  span  are  (utilizing  two-dimensional 
aerodynamic  oscillatory  coefficients) 


(39) 


'  1 

2  4 

r  " 

L‘  b 

IT  1)  U)  b  cos  .  J\. 

c/4  — 

h/b 

M’ 

►  --  < 

a 

T» 

,  j 

ft 

,  j 

where  the  matrix  -A.  contains  the  oscillatory  coefficients  as  given  in  Ref¬ 
erences  2  and  3,  that  is 


Hh 

^  a 

Lh/3 

Mah 

M 

O'O' 

Maf, 

V 

T 

lfta 

TIW 
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From  the  principle  of  virtual  work,  the  work  done  by  these  forces  acting 
through  a  virtual  displacement  is  equal  to, 


d  W 


/ 


L'  b 6—  +  M'  da  +  T'  d/3 
b 


dY 


(40) 


where  i  is  the  semi-span  parallel  to  the  Y-axis. 


Expressing  Equation  40  in  matrix  form, 


A. 

6W  =  /  |_  6b  6a  6/3 J 


L'b 

M' 

T* 


dY 


(41) 


Substituting  Equations  36,  37,  38,  and  39  into  Equation  41  results  in  a 
virtual  work  expression  in  terms  of  the  generalized  coordinates: 


4 

6W  =  63'  f  7rpo>2b4  cos*^^’  T(2)T(1)^dYq  (42) 


The  generalized  force  required  for  Equation  1,  Section  3. 1  is  simply, 


6W. 


QJ  = 


(43) 


lJ 


Applying  Equation  43  to  Equation  42  results  in  the  required  column  of 
generalized  aerodynamic  forces: 


Q  =  Trpco 


/  b4  cos  *c/4<P' '  T(1)*  T(2)'yv  T(2)  T{1)<p  dY  q 


or 


Q  =  cj  A  cj[ 


(44) 
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where  A  «  ». D  J  b4  cos  *c/<  £  T(1)'  T|2)'aT|Z)  T<L)  <p  dY 

o 

and  is  termed  the  generalized  aerodynamic  matrix. 

Aspect  Ratio  Corrections 


Aspect  ratio  corrections  were  applied  to  the  two-dimensional  aerody¬ 
namics  which  forced  the  air-load  distribution  to  take  on  an  elliptical 
distribution. 

The  two-dimensional  air  forces  used  were  seen  to  be  equal  to  the 
following: 


L'  b 

2  4 

np  U>  b  COS 

'  h/b ' 

M' 

.  =  < 

a 

T*  . 

k  P  . 

The  application  of  Equation  39  required  the  division  of  the  lifting  surface 
into  a  finite  number  of  strips,  evaluating  the  above  expression  at  n  sta¬ 
tions  and  integrating  numerically  across  the  span  to  obtain  the  generalized 
aerodynamic  matrix. 

With  the  above  representation  of  the  air  forces,  a  factor  was  determined 
which,  when  multiplied  by  Equation  39  at  each  station,  forced  the  air 
load  to  assume  the  required  elliptical  distribution. 

Consider  the  function  G  plotted  in  Figure  5,  in  which  the  function  assumes 
the  value  of  one  at  Y  =  0  and  zero  at  Y  =  l  and  describes  an  ellipse 
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L. 


The  equation  of  the  function  G  may  be  written  as 


or 


G  = 


‘•4 

i 1 


(45) 


Equation  45  lends  itself  to  a  factor  which  can  be  used  to  modify  the  air 
forces  given  in  Equation  39. 


Letting  G  =  F. 


-  JV- 


.2  ,.2 


Y  /f  ,  Equation  39  may  be  written  as 


L'b 

2  4 

7T  p  0)  b  cos  *  .  F„ 

c/4  F  — 

h/b 

M' 

►  =  i 

a 

,T’  . 

.  P  . 

(46) 


Following  the  samr  procedure  as  in  Paragraph  3.1.3,  the  column  of  gen¬ 
eralized  forces  is, 

q=irPo?  I  b4  cos  Ff  <£  T(1*'  T^T^^dYq 

o 


or 


2 

Q  =  u>  A  q 

where  A,  the  generalized  aerodynamic  matrix,  is  now  equal  to 

t 

n  p  f  b4  cos  *  F _<P'  T(1)'  dY 

J  c/4  F-  -  - - -  x 

o 


3.2 


SYMMETRIC  FLUTTER  ANALYSIS 


Study  of  the  flutter  characteristics  of  the  empennage,  symmetric  sense, 
was  initially  restricted  to  an  investigation  of  the  basic  horizontal  stabili¬ 
zer  with  an  assumed  cantilevered  restraint  at  the  horizontal  stabilizer 
pivot.  Subsequent  coupling  of  the  stabilizer  with  the  fuselage  allowed  an 
investigation  into  this  added  effect. 

3.2.1  Cantilevered  Analysis 

The  constrained  empennage  coordinates  are  the  appropriate  normal  vibra¬ 
tion  modes  of  the  horizontal  stabilizer  consistent  with  the  assumed  canti¬ 
levered  restraint.  Therefore,  the  following  generalized  coordinates  are 
considered: 


-  Horizontal  stabilizer  rigid  body  pitch  (about  pitch  axis) 

**2 

% 


-  Elevator  rotation 


-  1st  coupled  horizontal  stabilizer  elastic  mode 


q4  -  2nd  coupled  horizontal  stabilizer  elastic  mode 
q  -  3rd  coupled  horizontal  stabilizer  elastic  mode 


3 . 2 . 1 . 1  Generalized  Mass  Matrix 


Since  only  the  horizontal  tail  is  involved,  the  total  generalized  mass 
matrix  is: 


«  "  — H.T. 

For  the  generalized  coordinates  given  in  Paragraph  3.2.1,  the  generalized 
displacement  matrix  for  the  horizontal  tail  takes  the  following  form  for 
the  displacements  indicated  in  Equation  21,  Paragraph  3. 1. 1: 
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where  the  </>jj  are  column  matrices  relating  the  appropriate  column 
matrix  6  to  the  j1*1  generalized  coordinate. 

The  mass  matrix  J  for  the  horizontal  tail  is  that  given  by  Equation  21, 
Paragraph  3.1.1. 

3.2. 1.2  Aerodynamic  Matrix 

For  the  generalized  coordinates  stated  in  Paragraph  3.2.1,  the  general 
ized  displacement  matrix  ^is  that  given  in  Paragraph  3.2. 1. 1,  and  the 
generalized  uerodynamic  matrix  is  simply 


The  transformation  matrices  T  take  on  the  following  forms: 
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0  1.0  0  -a 

0  0  0  0  1.0 

0  0  0  0  0 


0 
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0 

1.0 


where  a  is  the  streamwise  distance  from  the  quarter  chord  to  the  elastic 
axis,  positive  aft  of  the  quarter  chord. 


3.2.2  Free-Free  Analysis 


The  unrestrained  fuselage  coordinates  considered  for  this  analysis  are 
the  free-free  normal  vibration  modes  of  the  fuselage  with  wing  and 
empennage  mass  effects  included.  The  empennage  coordinates  are  the 
same  normal  vibration  modes  utilized  in  the  cantilevered  analysis. 
Therefore,  the  following  generalized  coordinates  were  considered: 


% 


‘6 


m8 

% 


Aircraft  vertical  translation 
Aircraft  pitch 

Horizontal  stabilizer  rigid  body  pitch  (about  pitch  axis) 

Elevator  rotation 

1st  fuselage  vertical  bending  mode 

2nd  fuselage  vertical  bending  mode 

3rd  fuselage  vertical  bending  mode 

4th  fuselage  vertical  bending  mode 

1st  coupled  horizontal  stabilizer  elastic  mode 

2nd  coupled  horizontal  stabilizer  elastic  mode 


3.2.2. 1  Generalized  Mass  Matrix 

For  the  free-free  analysis,  the  total  generalized  mass  matrix  is  of  the 
form, 


M  =  2  M„  _  +  M_ 
—  -H.T.  — F 
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where  the  wing  mass  and  mass-inertia  effects  are  included  as  lumped 
elements  in  the  fuselage  mass  distribution. 

For  the  generalized  coordinates  given  in  Paragraph  3.2.2,  the  general¬ 
ized  displacement  matrix  (p  for  the  horizontal  tail  takes  the  following 
form  for  the  displacements  indicated  in  Equation  21,  Paragraph  3.1.1: 
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where  the  <p^  are  column  matrices  relating  the  appropriate  column 
matrixes  to  the  generalized  coordinate. 

The  fuselage  generalized  displacement  matrix  £  may  be  formulated  from 

Equation  14,  Paragraph  3.1.1,  by  considering  only  two  displacements, 

6  and  6  .  Therefore, 
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where  the  <p^  elements  are  again  column  matrices  relating  the  specified 
displacements  and  6^  to  the  stated  generalized  coordinates. 

The  mass  matrix  J  for  the  horizontal  tail  is  that  given  by  Equation  21, 
Paragraph  3.1.1. 
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For  the  fuselage,  since  only  two  displacements  are  considered,  the  mass 
matrix  J  takes  the  following  form: 


where  the  elements  also  reflect  the  vertical  tail  mass  and  mass-inertia 
effects . 

3. 2. 2. 2  Generalized  Aerodynamic  Matrix 

The  aerodynamic  matrix  suitable  for  this  analysis  is  simply 


where  the  transformation  matrices  T  are  as  given  in  Paragraph  3.2. 1.2 
and  the  generalized  displacement  matrix  <p  is  shown  in  Paragraph 
3.2.2. 1. 


3.3  ANTISYMMETRIC  FLUTTER  ANALYSIS 

Study  of  the  flutter  characteristics  of  the  empennage,  antisymmetric 
sense,  as  in  the  symmetric  case,  was  initially  restricted  to  an  investi¬ 
gation  of  the  basic  vertical  stabilizer  with  an  assumed  cantilevered 
restraint  at  the  vertical  stabilizer-fuselage  intersection.  This  initial 
analysis  also  considered  a  flexible  horizontal  stabilizer  with  an  assumed 
cantilevered  restraint  at  the  horizontal  stabilizer  pivot.  Subsequent 
coupling  of  the  empennage  with  the  fuselage  provided  a  set  of  new  coordi¬ 
nates  which  then  were  utilized  in  the  flutter  analysis. 

3.3.1  Cantilevered  Ana  lysis 

The  constrained  empennage  coordinates  are  the  appropriate  normal 
vibration  modes  of  the  horizontal  and  vertical  stabilizers  consistent  with 
the  assumed  cantilevered  restraint  as  mentioned  in  Paragraph  3.3. 
Therefore  the  following  generalized  coordinates  were  considered: 

-  Rudder  rotation 
<*2  - 


Elevator  rotation 


q^  -  Horizontal  stabilizer  rigid  body  flapping 

-  Horizontal  stabilizer  rigid  body  yawing 

q  -  1st  coupled  horizontal  stabilizer  elastic  mode 
5 

q  -  2nd  coupled  horizontal  stabilizer  elastic  mode 
q  -  3rd  coupled  horizontal  stabilizer  elastic  mode 


-  1st  coupled  vertical  stabilizer  elastic  mode 

% 

q 


-  2nd  coupled  vertical  stabilizer  elastic  mode 


10 


-  3rd  coupled  vertical  stabilizer  elastic  mode 


3.3. 1.1  Generalized  Mass  Matrix 

For  this  phase  of  the  antisymmetric  analysis,  the  generalized  mass 
matrix  is  of  the  form, 


—  2—  H.T.  +  —  V.T. 


For  the  generalized  coordinates  given  in  Paragraph  3.3.1,  the  general¬ 
ized  displacement  matrix  <p{or  the  horizontal  tail  takes  the  following 
form  for  the  displacements  indicated  in  Equation  21,  Paragraph  3.1.1. 
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The  generalized  displacement  matrix  for  the  vertical  tail  may  be 
written  in  the  following  form  for  the  same  displacements  as  forthe  hori¬ 
zontal  tail  except  for  the  displacement 6  which  is  replaced  by  6  : 
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The  mass  matrix  J  for  the  horizontal  tail  is  that  given  by  Equation  21, 
Paragraph  3.1.1. 

The  mass  matrix  J  for  the  vertical  tail  is  given  by  Equation  22,  Para¬ 
graph  3.1.1  with  the  subscript  Z  replaced  by  Y. 


3.3. 1.2  Generalized  Aerodynamic  Matrix 

Since  we  have  aerodynamics  acting  on  both  the  horizontal  and  vertical 
tails,  the  total  aerodynamic  matrix  is 


A 


+  A 


V.T. 


where  the  transformation  matrices  for  both  the  horizontal  and  vertical 
tails  are  as  given  in  Paragraph  3.2. 1.2. 

The  generalized  displacement  matrices  ^  for  both  the  horizontal  and 
vertical  tails  are  as  given  in  Paragraph  3.3. 1. 1. 

3.3.2  Free- Free  Analysis 

Due  to  the  large  number  of  coordinates  required  to  describe  adequately 
the  motion  of  the  aircraft  in  an  antisymmetric  sense,  a  set  of  new  gen¬ 
eralized  coordinates  was  determined  for  use  in  the  flutter  analysis. 


3.3.2. 1  Modal  Coupling 


For  determination  of  a  new  set  of  coordinates,  the  unrestrained  fuselage 
coordinates  considered  are  the  free-free  normal  vibration  modes  of  the 
fuselage  with  wing  and  empennage  mass  and  mass-inertia  effects  included. 
The  empennage  coordinates  are  the  same  normal  vibration  modes  utilized 
in  the  cantilevered  analysis.  The  following  generalized  coordinates  were 
considered  for  the  modal  coupling  phase: 

q^  -  Aircraft  side  translation 

q  -  Aircraft  yaw 

q  -  Aircraft  roll 

u 

q^  -  1st  fuselage  lateral  bending  mode 

q5  -  2nd  fuselage  lateral  bending  mode 

q  -  3rd  fuselage  lateral  bending  mode 

q7  -  4th  fuselage  lateral  bending  mode 

q  -  1st  fuselage  torsion  mode 
8 

-  2nd  fuselage  torsion  mode 

q,rt  -  3rd  fuselage  torsion  mode 
10 

^11  ”  ^  fuselage  torsion  mode 
q  -  1st  coupled  vertical  stabilizer  elastic  mode 

1  u 

q  -  2nd  coupled  vertical  stabilizer  elastic  mode 

lu 

^14  ”  coupled  vertical  stabilizer  elastic  mode 

q  -  1st  coupled  horizontal  stabilizer  elastic  mode 

15 

q  -  2nd  coupled  horizontal  stabilizer  elastic  mode 

16 

^17  ~  3rt*  couPlet*  horizontal  stabilizer  elastic  mode 
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Rudder  rotation 

* 

Elevator  rotation 

Horizontal  stabilizer  rigid  body  flapping 

The  applicable  equations  of  motion  are  those  given  in  Equation  6,  Section 
3.1,  with  the  generalized  aerodynamic  matrix  A  equal  to  zero  and  the 
complex  eigenvalue  ft  reducing  to  G  =  1/^2  for  free  undamped  vibrations. 

With  the  elements  of  the  generalized  stiffness  matrix  K  taken  as  in  Equa¬ 
tion  35,  Paragraph  3. 1.2,  only  inertia  coupling  is  involved  in  the  genera¬ 
tion  of  the  new  set  of  generalized  coordinates. 

Generalized  Mass  Matrix 


Since  we  are  dealing  with  three  components,  i.e. ,  horizontal  and  verti¬ 
cal  tails  plus  the  fuselage,  the  generalized  mass  matrix  may  be  broken 
up  as  follows: 

M  =  2MH  T  +Mv  t  +Mf 

The  generalized  displacement  matrix  <p  for  the  horizontal  tail  for  the 


coordinates  given  in 

Paragraph  3.3.2. 

1  may  be  written  as 
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for  the  displacements  given  in  Paragraph  3. 1.1. 


For  the  vertical  tail,  the  generalized  displacement  matrix  </>  is 
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for  the  same  displacements  as  for  the  horizontal  tail  except  for  the  dis¬ 
placement  6  ,  which  is  replaced  by  6 
z  * 


The  generalized  matrix  <p  for  the  fuselage  is 

*11  *12  *13  *14  *15  *10  *17  °  °  °  °  0  0  0  0  0  0  0  0  0 

°  °  *23  °  °  °  °  *28  *20  *210  *211  0  0  0  0  0  0  0  0  0 

°  *32  0  *34  *35  *30  *37  °  °  °  0  0  0  0  0  0  0  0  0  0 


for  only  thrr  <  nsidered  displacements,  6 


<5  and 
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The  mass  matrix  J_  for  the  horizontal  tail  is  that  given  by  Equation  21, 
Paragraph  3.1.1. 

The  mass  matrix  J  for  the  vertical  tail  is  given  by  Equation  22,  Para¬ 
graph  3.1.1  with  the  subscript  Z  replaced  by  Y. 


For  the  three  displacements  considered  for  the  fuselage,  the  mass 
matrix  J  takes  the  following  form: 

r »  -sx  oi 


L°  0  y 

where  the  elements  also  reflect  the  wing  mass  and  mass-inertia  effects. 


3. 3. 2. 2  Flutter  Analysis 

The  new  set  of  generalized  coordinates  determined  from  a  solution  of  the 
equations  of  motion  described  in  Paragraph  3.3.2. 1,  is  as  follows: 

q^  -  Aircraft  side  translation 

-  Aircraft  yaw 

qg  -  Aircraft  roll 

q^  -  Horizontal  stabilizer  rigid  body  yawing 

q  -  1st  coupled  fuselage  -  empennage  elastic  mode 

5 

q  -  2nd  coupled  fuselage  -  empennage  elabtic  mode 

q^  -  3rd  coupled  fuselage  -  empennage  elastic  mode 

q  -  4th  coupled  fuselage  -  empennage  elastic  mode 

O 

q^  -  5th  coupled  fuselage  -  empennage  elastic  mode 

q^0  -  6th  coupled  fuselage  -  empennage  elastic  mode 

Generalized  Mass  Matrix 

As  in  the  modal  coupling  phase,  the  generalized  mass  matrix  may  be 
split  up  into  parts: 

2mh.t.  +  **v.t.  +  «F 

where  the  generalized  displacement  matrix  <p  for  the  horizontal  tail, 
vertical  tail  and  fuselage  are,  respectively, 
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The  mass  matrices  J.  for  the  three  components  are  as  given  in  Paragraph 
3.3.2. 1. 
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Generalized  Aerodynamic  Matrix 


As  in  the  antisymmetric  cantilevered  analysis,  the  aerodynamic  matrix 
may  be  written  as 


A  = 


+  A 


V.T. 


where  the  transformation  matrices  for  both  the  horizontal  and  vertical 
tail  are  as  given  in  Paragraph  3.2. 1.2. 

The  generalized  displacement  matrix  <p  for  both  the  horizontal  and  verti 
cal  tail  are  given  under  Generalized  Mass  Matrix. 


4.0  DISCUSSION  AND  RESULTS 


The  presentation  and  discussion  of  the  results  of  the  empennage  flutter 
analysis,  in  general,  follows  the  sequence  of  Sections  3.2  and  3.3. 

Also,  this  section  presents  the  idealization  of  the  empennage  used  for  the 
analysis,  the  structural  properties  of  such  structure,  and  other  data 
required  for  the  flutter  analysis  of  the  empennage. 

4.1  IDEALIZATION  OF  AIRCRAFT 

As  previously  mentioned,  structural  idealization  of  the  empennage  and 
fuselage  was  restricted  to  equivalent  beams  possessing  both  flexure  and 
torsional  stiffness  distributions,  with  the  aerodynamic  loading  repre¬ 
sented  by  two-dimensional  or  strip-theory  aerodynamic  forces. 

Figure  6  of  Reference  1  presents  the  actual  structure  of  the  XV-5A 
empennage  and  fuselage  to  be  idealized. 

4.1.1  Structural  Representation 

As  the  aircraft  structure  shown  in  Figure  6  of  Reference  1  shows,  this 
may  be  thought  of  as  a  number  of  components  assembled  to  form  the 
complete  structure,  i.e.,  fuselage  and  horizontal  and  vertical  stabili¬ 
zers  plus  control  surfaces.  From  this  standpoint,  the  following  para¬ 
graphs  deal  with  these  components  as  they  were  treated  during  the 
analysis. 

4. 1.1.1  Horizontal  Stabilizer- Elevator 

Figure  6  shows  the  planform  of  the  horizontal  stabilizer  and  elevator. 

The  elastic-axis  shown  represents  the  equivalent  beam  which  replaces 
the  actual  structure  shown  in  Figure  6  of  Reference  1.  The  elevator 
was  taken  as  a  rigid  body  having  distributed  mass  and  mass-inertia 
properties.  The  elevator  system  is  manual;  therefore,  its  rotational 
restraint  was  idealized  into  an  equivalent  spring.  This,  in  conjunction 
with  the  pitching  inertia  of  the  elevator,  resulted  in  an  elevator  rota¬ 
tional  frequency,  which  was  varied  to  cover  expected  ranges.  Figure 
8  presents  a  plot  of  the  spanwise  distribution  of  the  bending  and 
torsional  stiffnesses  for  the  horizontal  stabilizer.  These  stiffness 
values  reflect  a  maneuver  flight  condition,  i.e.,  they  reflect  a  mini¬ 
mum  stiffness  due  to  skin  ineffectiveness. 
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4 . 1 . 1 . 2  Vertical  Stabilize r- Rudder 


Figure  7  shows  the  planform  of  the  vertical  stabilizer  and  rudder,  with 
the  elastic-axis  shown  representing  the  orientation  of  the  equivalent 
beam  which  replaces  the  actual  structure  shown  in  Figure  6  of 
Reference  1.  The  rudder  system,  being  equivalent  to  the  elevator 
system,  was  Idealized  in  the  same  fashion.  Figure  9  presents  a  plot 
of  the  spanwise  distributions  of  the  bending  and  torsional  stiffnesses 
for  the  vertical  stabilizer.  As  for  the  horizontal  stabilizer,  these 
plots  reflect  a  minimum  stiffness,  also. 

4. 1.1.3  Fuselage 

The  fuselage  structure  shown  in  Figure  6  of  Reference  1  was  idealized 
into  an  equivalent  beam  having  vertical,  lateral  and  torsional  stiffness. 
These  distributions  are  shown  in  Figure  10.  The  elastic  axis  of  the 
fuselage  was  assumed  to  be  at  the  intersection  of  the  vertical  plane  at 
B.  L.  O  with  the  horizontal  plane  at  W.  L.  113.00.  With  this  represen¬ 
tation,  the  vertical  stabilizer  elastic-axis  intersects  the  fuselage  at 
F.S.  459.37,  as  shown  in  Figure  7. 

4.1.2  Aerodynamic  Representation 

Since  aerodynamic  forces  acting  on  the  wing  normally  are  small,  com¬ 
pared  to  forces  on  the  empennage  when  predominant  empennage  motion 
is  involved,  the  wing  contribution  to  the  formulation  of  the  aerodynamic 
matrix  was  neglected.  Flutter  analyses  were  performed  for  several 
different  altitudes  to  observe  effects  of  air  density.  The  aerodynamic 
coefficients  used  reflected  both  incompressible  and  compressible  flow. 
The  former  coefficients  were  based  largely  on  the  work  of  References  2 
and  3,  whereas  the  compressible  coefficients  were  derived  from  those 
covered  by  References  4  and  5. 

4. 1.2.1  Empennage 

Application  of  Equation  44,  Paragraph  3. 1.3,  required  the  division  of 
the  horizontal  and  vertical  tails  into  a  number  of  spanwise  strips 
oriented  parallel  to  the  free  stream.  Evaluation  of  Equation  44  at  each 
of  the  air-force  stations  was  then  required,  and  the  total  aerodynamic 
matrix  for  the  surface  was  formed  by  numerical  integration.  Figures  6 
and  7  also  depict  the  planforms  which  were  used  for  generation  of  the 
aerodynamic  matrix  for  the  horizontal  and  vertical  tails,  respectively. 


4.1.3 


Mass  Representation 


Only  one  gross  weight  condition,  that  of  the  design  gross  weight  of 
approximately  9200  lbs,  was  studied  during  the  analysis.  This  weight 
condition  was  applied  only  to  the  free-free  analysis,  both  symmetric 
and  antisymmetric.  From  a  given  mass  distribution,  the  distribution 
was  reworked  to  reflect  a  lumped-parameter  system  which  was 
governed  by  an  axis  system  dictated  by  Equation  21  or  22  of  Paragraph 
3.1.1.  Inconsistencies  between  distributions  shown  for  the  symmetric 
and  antisymmetric  analyses  are  attributed  to  the  relative  time  between 
completion  of  each  analysis  with  the  consequent  up-dating  of  a  distri¬ 
bution. 

4 . 1 . 3 . 1  Symmetric  Analysis 

The  mass  distribution  of  the  two  components,  the  horizontal  tail  and 
fuselage,  required  for  both  the  initial  cantilever  analysis  and  the  free- 
free  analysis,  are  presented  in  Tables  1  and 2  for  the  horizontal  tail  and 
fuselage,  respectively.  The  properties  shown  reflect  the  nomenclature 
as  described  by  Equation  21  of  Paragraph  3.1.1.  Elevator  and  elevator 
mass-balance  are  also  shown  separately  to  indicate  the  mass-balance 
distribution  utilized  in  the  analysis.  Table  2  presents,  in  addition  to 
the  basic  fuselage  distribution,  the  wing  and  empennage  mass  proper¬ 
ties  which  made  up  the  inertial  force  and  moment  contributions  of 
these  components  to  the  fuselage  symmetric  normal  modes. 

4. 1.3.2  Antisymmetric  Analysis 


The  mass  distribution  of  the  three  components,  namely  the  horizontal 
tail,  vertical  tail  and  fuselage,  required  for  both  the  initial  cantilevered 
analysis  and  the  free-free  analysis,  are  presented  in  Tables  3  through 
5  for  the  horizontal  tail,  vertical  tail  and  fuselage,  respectively.  The 
properties  shown  again  reflect  the  nomenclature  as  described  by  Equa¬ 
tion  21  of  Paragr  aph  3.1.1.  To  show  elevator  and  rudder  mass-balance, 
separate  tables  of  the  balanced  control  surface  and  of  the  balance  weights 
are  also  shown  in  Tables  3  and  4.  Also  shown  in  Table  4  are  the  hori¬ 
zontal  tail  mass  properties  which  made  up  the  inertial  force  and  moment 
contribution  of  this  component  to  the  vertical  stabilizer  normal  modes. 
Table  5  presents,  in  addition  to  the  basic  fuselage  distribution,  the 
wing  and  empennage  mass  properties  which  made  up  the  inertial 
force  and  moment  contributions  of  these  components  to  the  fuselage 
antisymmetric  normal  modes. 
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4.2 


MODAL  FUNCTIONS 


As  mentioned  in  Section  3.0,  the  equations  of  motion  were  formulated 
in  terms  of  modal  functions.  These  were  principally  the  normal  modes 
of  vibration  of  the  component  consistent  with  an  assumed  constraint. 

4.2.1  Symmetric  Analysis 

As  mentioned  in  Section  3.2,  the  components  involved  in  the  symmetric 
analysis  were  the  horizontal  stabilizer,  fuselage  and  elevator.  The 
horizontal  stabilizer  normal  modes  are  the  coupled  modes  (inertial)  of 
bending  and  torsion  of  the  equivalent  beam  considered  cantilevered  at 
the  horizontal  stabilizer  pivot.  The  normal  modes  of  the  fuselage  are 
the  uncoupled  free-free  modes  of  the  fuselage  equivalent  beam  with 
wing  and  empennage  mass  and  mass-inertia  effects  included.  Vertical 
shear  effects  (AG)  are  also  considered  in  the  determination  of  these 
latter  modes. 

4.2. 1.1  Cantilevered 

Figures  11,  12  and  13  depict  the  first  three  cantilevered  mode  shapes 
of  the  horizontal  stabilizer,  white  Table  6  gives  a  tabulation  of  the 
modal  data.  The  vertical  bending  is  denoted  by  (/’whereas  the  torsion 
of  the  surface  about  the  elastic-axis  is  denoted  by  0. 

4. 2. 1.2  Free-Free 


In  addition  to  the  horizontal  stabilizer  normal  modes,  the  fuselage 
degrees  of  freedom  were  included,  and  these  shapes  are  plotted  in 
Figures  14  and  15,  with  the  modal  data  tabulated  in  Table  7. 

4.2.2  Antisymmetric  Analysis 

Also,  as  mentioned  in  Section  3.3,  the  components  involved  in  the  anti¬ 
symmetric  analysis  were  the  horizontal  stabilizer,  vertical  stabilizer, 
and  fuselage,  plus  elevator  and  rudder.  The  horizontal  stabilizer  nor¬ 
mal  modes  are  the  coupled  modes  (inertial)  of  bending  and  torsion  of 
the  equivalent  beam  considered  cantilevered  at  the  horizontal  stabilizer 
pivot.  The  vertical  stabilizer  normal  modes  are  also  the  inertially 
coupled  modes  of  bending  and  torsion  of  the  equivalent  beam  considered 
cantilevered  at  the  intersection  of  the  elastic-axis  with  the  assumed 
fuselage  elastic-axis  (F.S.  459.37).  Mass  and  mass-inertia  effects 
of  the  horizontal  tail  were  considered  in  these  latter  modes.  The 


normal  modes  of  the  fuselage  are  the  uncoupled  free-free  modes  of  the 
fuselage  equivalent  beam,  with  wing  and  empennage  mass  and  mass- 
inertia  effects  included. 

4.2.2. 1  Cantilevered 

Figures  16,  17  and  18  present  the  first  three  cantilevered  mode  shapes 
of  the  horizontal  stabilizer,  while  Table  8  tabulates  the  modal  data. 

The  equivalent  information  for  the  vertical  stabilizer  is  presented  in 
Figures  19,  20,  and  21  and  Table  9. 

4. 2. 2. 2  Free-Free 

For  the  free-free  analysis,  as  mentioned  in  Paragraph  3.3.2,  coupled 
modes  for  the  aircraft  were  calculated,  and  a  new  set  of  generalized 
coordinates  was  chosen  for  the  flutter  analysis.  The  horizontal  and 
vertical  stabilizer  normal  modes  utilized  in  the  coupling  are  the  same 
constrained  modes  of  vibration  utilized  in  the  cantilevered  analysis. 

The  fuselage  mode  shapes  are  presented  in  Figures  22  and  23  for  the 
uncoupled  lateral  bending  and  in  Figures  24  and  25  for  uncoupled 
torsion,  the  modal  data  being  shown  in  Table  10. 

4,3  FLUTTER  ANALYSIS 

As  explained  previously,  the  analysis  of  the  empennage  has  been 
divided  into  a  symmetric  and  an  antisymmetric  investigation.  Each  of 
these  analyses  was  then  divided  further  into  basic  cantilevered  analyses 
and  then  expended  to  allow  for  fuselage  flexibility  and  aircraft  rigid 
body  degrees  of  freedom . 

In  general,  most  of  the  results  of  the  flutter  analyses  are  plotted  in 
the  conventional  sense,  i.e.,  overall  damping  (g)  versus  forward 
velocity  (V).  In  addition,  the  flutter  frequency  (f)  is  also  plotted 
versus  velocity.  Equation  6  of  Section  3 . 1  implies  that  these  will  be 
n-r  roots,  where  n  is  the  number  of  generalized  coordinates  considered 
for  the  problem,  and  r  is  the  number  of  generalized  coordinates  de¬ 
fining  rigid  body  degrees  of  freedom  of  zero  frequency.  Since  genera¬ 
tion  of  the  air-force  matrix  A  requires  the  selection  of  a  reduced  fre¬ 
quency  parameter  (kQ),  the  curves  of  damping-frequency  versus  flutter 
speed  were  generated  by  considering  a  range  of  reduced  frequency 
parameters.  The  flutter  parameters  were  noted  by  recording  the  lowest 
velocity  at  which  a  V-g  plot  crossed  the  g  =  o  axis  and  then  noting  the 
flutter  frequency  at  this  velocity  from  the  V-f  plot  for  the  appropriate 
root. 
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4.3.1  Symmetric 


The  initial  phase  of  the  symmetric  analysis,  as  explained  previously, 
was  restricted  to  the  basic  horizontal  tail,  which  was  assumed  canti¬ 
levered  at  the  pivot.  The  next  phase  included  fuselage  effects  and  also 
aircraft  rigid  body  degrees  of  freedom. 

4.3. 1.1  Cantilevered  Analysis 

The  large  number  of  parameters  involved  in  the  flutter  analysis,  such 
as  elevator  rotational  frequency,  horizontal  tail  pitch  frequency,  density, 
aerodynamic  parameters  (aspect-ratio  cor  rections  and  effect  of  aero¬ 
dynamic  balance)  precluded  an  extensive  variation  of  these  parameters. 
Therefore,  only  those  parameters  which  were  felt  to  be  significant 
were  varied  over  a  larger  range .  Both  incompressible  and  compressible 
aerodynamic  coefficients  were  used  in  this  analysis,  as  the  incompress¬ 
ible  analysis  indicated  that  a  flutter  problem  existed  on  the  basic 
horizontal  tail. 

Incompressible 


Figures  26  through  30  present  the  results  of  the  basic  cantilevered 
analysis  for  a  variable  horizontal  tail  pitch  frequency  of  25  to  80  cycles 
per  second,  for  a  constant  elevator  rotational  frequency  of  10  cycles 
per  second  and  at  an  altitude  of  sea-level.  Aerodynamic  parameters 
were  such  that  aspect-ratio  corrections  and  effects  of  aerodynamic 
balance  were  neglected. 

Figures  31  through  33  present  results  equivalent  to  those  in  Figures  26 
through  30,  but  for  an  elevator  rotational  frequency  of  30  cps  and  for 
horizontal  tail  pitch  frequencies  of  40,  60  and  80  cps. 

Figures  34  through  36  show  the  results  of  removing  the  elevator  degree 
of  freedom  (assuming  elevator  to  be  rigidly  clamped)  from  the  analysis. 
Other  parameters  are  as  in  Figures  31  through  33. 

Figures  37  through  39  present  effects  of  aspect-ratio  corrections 
(tip- round -off)  for  variable  horizontal  tail  pitch  frequencies  of  40,  60 
and  80  cps,  for  an  elevator  rotational  frequency  of  10  cps,  and  at 
sea- level  altitude. 

Figure  40  shows  the  effect  of  including  aerodynamic  balance  in  the  air¬ 
force  matrix  for  a  horizontal  tail  pitch  frequency  of  60  cps  and  an 
elevator  rotational  frequency  of  10  cps,  all  at  sea-level  altitude. 


Aerodynamic  balance  effects  were  accounted  for  by  means  of  the  methods 
described  in  References  2  and  3.  These  references  treat  any  balance 
as  external,  (nose  overhang  exposed  directly  to  the  free  stream), 
whereas  the  actual  aerodynamic  balance  in  the  elevator  is  of  an  internal 
nature.  It  was  felt  that  treatment  of  the  balance  effects  by  means  of 
the  techniques  of  References  2  and  3  would  indicate  trends,  which  would 
be  sufficient  for  the  analysis. 

Figures  41  through  43  present  altitude  effects  for  a  pitch  frequency  of 
60  cps  and  an  elevator  rotational  frequency  of  10  cps.  An  altitude  of 
9300'  represents  the  minimum  altitude  at  which  the  maximum-limit 
dive  Mach  number  is  reached. 

A  cross-plot  of  horizontal  tail  pitch  frequency  versus  flutter  speed  is 
presented  in  Figure  44  for  the  different  parameteric  studies  shown  in  the 
V-g  plots  of  Figures  26  through  43. 

Compressible 

Compressible  flutter  analyses,  using  aerodynamic  coefficients  for  a 
Mach  number  of  0.60,  0.70  and  0.80,  were  conducted  for  a  fixed 
horizontal  tail  pitch  frequency  of  60  cps  and  an  elevator  rotational  fre¬ 
quency  of  10  cps.  Aerodynamic  parameters  (aspect-ratio  corrections 
and  effects  of  aerodynamic  balance)  were  neglected  for  this  phase  of 
the  analysis,  with  altitude  variations  taken  for  the  four  altitudes  used 
in  the  incompressible  analysis. 

Figures  46  through  48  present  the  results  for  the  sea-level  condition, 
while  the  results  for  the  9300*,  20,  000'  and  30,000'  analyses  are  pre¬ 
sented  in  Figures  49  through  57.  Only  the  roots  of  concern  have 
been  plotted  on  both  the  v-g  and  V-f  plots. 

The  flutter  speed  deduced  from  the  above  figures  represents  an 
apparent  speed  and  not  a  true  flutter  speed.  This  is  because  the  speed 
corresponding  to  the  particular  Mach  number,  in  general,  is  not 
equal  to  the  apparent  flutter  speed.  For  this  reason,  an  auxiliary  plot 
is  given  in  Figure  58,  wherein  the  apparent  flutter  speeds  deduced 
from  Figures  46  through  57  have  been  converted  to  equivalent  air  speeds. 
Crossing  of  a  particular  altitude  line  with  the  plot  generated  from  the 
apparent  flutter  speeds  then  yields  a  true  flutter  speed.  Also  included 
in  Figure  58  are  the  incompressible  results  obtained  from  Figures  29 
and  41  through  43. 
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4.3. 1.2  Free-Free  Analysis 


The  free-free  analysis  was  conducted  for  two  values  of  the  elevator 
rotational  frequency  and  for  three  values  of  the  horizontal  tail  pitch 
frequency.  The  analysis  was  conducted  for  a  sea-level  condition  with 
the  aerodynamic  parameter  (aspect-ratio  corrections  and  effects  of 
aerodynamic  balance)  neglected.  Figures  59  through  61  present  the 
results  of  the  analysis  for  an  elevator  rotational  frequency  of  10  cps, 
while  the  same  equivalent  results  for  an  elevator  rotational  fre¬ 
quency  of  infinity  (locked  elevator)  are  presented  in  Figures  62  through 
64,  with  a  cross  plot  of  the  overall  effects  shown  in  Figure  65. 

4.3.2  Antisymmetric 


The  initial  phase  of  the  antisymmetric  analysis  was  restricted  to  the 
basic  vertical  tail,  assumed  cantilevered  at  the  fuselage-vertical 
stabilizer  intersection.  This  analysis  included  the  effects  of  a  flex¬ 
ible  horizontal  stabilizer.  The  second  phase  included  fuselage  effects 
and  also  aircraft  rigid  body  degrees  of  freedom. 

4.3.2. 1  Cantilevered  Analysis 

As  in  the  symmetric  analysis,  the  large  number  of  parameters  in¬ 
cluded  in  the  flutter  analysis  of  the  empennage,  such  as  elevator  and 
rudder  rotational  frequency,  horizontal  tail  yawing  and  rolling,  density, 
and  aerodynamic  parameters  (aspect-ratio  corrections  and  effect  of 
aerodynamic  balance)  precluded  an  extensive  variation  of  these  para¬ 
meters.  Therefore,  only  the  parameters  considered  significant  were 
varied.  Only  incompressible  aerodynamic  coefficients  were  used  in  this 
analysis,  with  aspect  ratio  effects  and  aerodynamic  balancing  effects 
ignored.  The  end-plating  of  the  vertical  tall  by  the  tip-mounted,  hori¬ 
zontal  tail  fully  justified  the  use  of  the  two-dimensional  strip  theory,  while 
the  neglect  of  the  aerodynamic  balance  was  felt  to  be  of  small  concern  in 
the  overall  analysis. 

Figures  66  and  67  present  the  results  of  the  cantilevered  analysis  of 
the  empennage  for  a  rigid  body  yawing  frequency  of  10  and  20  cps. 

As  explained  in  Paragraph  4.3.2,  the  effects  of  a  flexible  horizontal 
stabilizer  were  included,  in  addition  to  the  rigid  body  modes  of  flapping 
(rocking)  and  yawing.  Antisymmetric  motion  of  the  elevators  was 
included,  with  the  frequency  taken  as  15  cps.  Variations  of  the  fre¬ 
quency  of  the  latter  mode,  as  shown  in  initial  analyses  not  presented 
in  this  report,  indicated  no  change  in  flutter  mode  or  speed.  Similar 
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earlier  studies  of  rudder  frequency  variations  showed  the  same  results; 
therefore,  the  rudder  rotational  frequency  was  held  to  10  cps  for  the 
analysis.  Horizontal  tail  flapping  was  held  to  10  cps,  which  was  felt 
to  be  reasonable  for  the  large  rolling  moment  of  inertia  involved. 

Studies  of  variations  of  this  parameter  in  initial  analyses  also  indicated 
no  change  in  flutter  speed  or  mode. 

4.3. 2. 2  Free-Free  Analysis 

As  explained  in  Paragraph  3.3.  2,  because  of  the  large  number  of  coor¬ 
dinates  required  to  define  adequately  the  antisymmetric  motion  of  the 
empennage,  which  would  include  effects  of  fuselage  stiffness  and  also 
aircraft  rigid  body  modes,  it  is  a  practical  impossibility  to  vary  all 
unknown  parameters.  Therefore,  certain  values  of  these  parametei'6 
were  fixed,  based  on  initial  analyses  of  the  empennage  assumed  canti¬ 
levered  from  the  fuselage.  With  this  in  mind,  several  analyses  were 
performed  for  two  values  of  the  horizontal  tail  yawing  frequency.  All 
other  structural  parameters,  such  as  rudder  and  elevator  rotational 
frequency  and  horizontal  tail  flapping,  were  held  to  certain  values  in  the 
modal  coupling  for  the  purpose  of  determining  the  set  of  new  generalized 
coordinates  required  for  the  flutter  analysis.  As  in  the  cantilevered 
analysis,  only  incompressible  aerodynamic  coefficients  were  used, 
with  aspect  ratio  effects  and  aerodynamic  balancing  effects  ignored. 

Figures  68  and  69  present  the  results  of  the  free-free  antisymmetric 
analysis  of  the  empennage  for  a  horizontal  tail  yawing  frequency  of  10 
and  20  cps,  respectively.  The  parameters  held  fixed  were  rudder 
rotation  of  10  cps,  elevator  rotation  (assumed  antisymmetric  motion) 
of  15  cps,  and  a  horizontal  tail  flapping  frequency  of  10  cps.  The 
latter  parameters  were  included  in  the  modal  coupling  and  were  con¬ 
sequently  reflected  in  the  six  antisymmetric  fuselage-empennage  modes 
used  as  generalized  coordinates  for  the  flutter  analysis. 
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5.0  CONCLUSIONS 


Because  of  the  manner  in  which  the  preliminary  flutter  analysis  of  the 
empennage  was  conducted,  i.e. ,  a  breaking  down  into  initial  component 
analysis,  conclusions  may  best  be  drawn  based  on  the  results  presented 
in  Section  4.0. 

5.1  SYMMETRIC  ANALYSIS  -  CANTILEVERED  CONDITION 

Figures  26  through  58  present  the  V-g  plots  of  these  analyses,  in 
addition  to  cross-plots  for  over  all  evaluation.  These  initial  analyses, 
preliminary  in  nature  and  done  at  a  time  when  only  parametric  studies 
can  be  made.  Indicate  (Figure  44)  that  the  controlling  parameter  on  the 
flutter  speed  of  the  horizontal  tail  is  that  of  the  horizontal  tail  rigid 
body  pitch  frequency. 

Figure  58  presents  the  proposed  flight  boundary  of  the  XV-5A  aircraft; 
a  15%  flutter  margin  on  the  limit  dive  speed  at  sea- level  is  575  knots. 
The  15%  requirement  on  therefore  necessitates  a  horizontal  tail 
pitch  frequency  of  approximately  55  cps,  based  on  incompressible 
theory.  The  trend  of  the  plots  shown  in  Figure  44  shows  that,  as  the 
pitch  frequency  is  increased  from  a  low  value  to  the  upper  value,  the 
flutter  speed  approaches  the  basic  bending  torsion  flutter  speed  of  the 
surface.  The  mode  of  flutter  involved  is  that  of  horizontal  tail  pitch 
coupled  with  the  first  coupled  elastic  mode  of  the  horizontal  tail.  The 
mode  shapes  of  this  mode  are  presented  in  Figure  11.  It  should  be 
noted  that  this  mode  exhibits  both  a  bending  and  torsion  component, 
the  shapes  of  which  are  characteristic  of  the  first  uncoupled  bending 
and  torsion  modes  of  a  cantilever  beam.  At  60  cps  pitch  frequency, 
the  flutter  frequency  is  40  cps.  Referring  to  Figure  29,  the  V-g  plot 
for  this  condition  shows  that  the  3rd  root  is  the  flutter  mode.  Tracing 
this  plot  back  to  zero  airspeed,  we  find  that  the  coupled  mode  of 
vibration  is  approximately  63  cps.  Solution  of  the  equations  of  motion 
for  zero  aerodynamics  yields  approximately  the  same  results  (65  cps). 
This  result  indicates  a  coupling  between  the  pitch  mode  and  the  torsion 
component  of  the  first  coupled  horizontal  tail  elastic  mode.  This  effect 
is  shown  in  Figure  45.  Note  that  the  frequencies  of  both  the  coupled 
and  uncoupled  elastic  modes  of  the  horizontal  stabilizer  are  plotted 
against  the  results  of  coupling  the  elastic  modes  with  a  variable  hori¬ 
zontal  tail  pitch  frequency,  for  a  constant  elevator  rotational  frequency 
of  10  cps.  From  this  plot,  it  can  be  seen  that  the  mode  of  flutter  of 
the  horizontal  tail  would  be  dependent  on  the  available  horizontal  tail 


pitch  frequency.  For  a  low  value  of  pitch  frequency  the  bending  mode 
would  predominate,  whereas,  for  a  higher  value  of  pitch  frequency,  the 
torsion  mode  would  be  predominant.  This  effect  is  evident  also  in 
Figures  26  through  39 . 

Variations  of  the  elevator  rotational  frequency  of  10  cps,  30  cps,  and 
finally  locking  out  of  the  elevator  (infinite  restraint)  showed  no 
appreciable  change  in  the  flutter  trends  of  the  surface. 

Changes  in  the  aerodynamic  parameters  shown  in  Figure  44  also  showed 
no  appreciable  change  in  the  basic  flutter  trend. 

Altitude  trends  are  shown  mainly  to  present  the  data.  These  indicate 
increasing  flutter  speeds,  as  is  normal. 

Compressible  analyses  were  conducted,  since  the  design  flight  envelope 
(see  Figure  58)  encompasses  the  high  subsonic  region  and  the  initial 
incompressible  analyses  indicated  that  a  relatively  high-pitch  frequency 
was  required.  Figure  58  presents  the  results  of  these  analyses  for  an 
assumed  horizontal  tail  pitch  frequency  of  60  cps.  The  M  =  o  results 
from  Figure  29  and  Figures  41  through  43  are  also  plotted.  Tlv'  trend 
of  the  plots  indicates  the  usual  drop  in  flutter  speed  as  the  transonic 
region  is  approached.  Crossing  of  a  flutter  plot  for  a  particular 
altitude  with  the  altitude  line  radiating  from  the  origin  yields  only  one 
or  two  true  flutter  speeds  (extrapolated  for  the  sea-level  and  9300' 
altitude).  The  boundary  defined  by  these  points  lies  within  the  15% 
margin  on  equivalent  airspeed,  which  indicates  that  a  possible  flutter 
situation  might  exist  on  the  XV-5A  aircraft  if  the  design  flight 
envelope  is  to  be  met. 

5 . 2  SYMMETRIC  ANALYSIS  -  FREE-FREE  CONDITION 

Figures  59  through  65  present  the  V-g  plots  of  these  analyses,  along 
with  a  cross-plot  for  various  elevator  rotational  frequencies,  for  three 
values  of  the  horizontal  tail  pitch  frequency.  The  lower  boundary  of 
Figure  65  (for  an  elevator  frequency  of  10  cps)  is  relatively  insensitive 
to  horizontal  tail  pitch  frequency,  whereas  the  upper  boundary  exhibits 
the  same  trend  as  for  the  cantilevered  analysis  (Figure  44).  The  flutter 
mode  for  the  lower  boundary  is  fuselage  bending  coupled  with  elevator 
rotation.  However,  raising  the  damping  level  by  0.03  (from  g  =  o), 
eliminates  this  boundary,  and  the  flutter  mode  is  now  in  effect, 
essentially  the  same  as  shown  for  the  cantilevered  case.  That  is,  the 
flutter  mode  is  horizontal  tail  pitch  coupled  with  predominant  horizontal 
tail  torsion  for  the  higher  pitch  frequencies. 
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5.3 


ANTISYMMETRIC  ANALYSIS  -  CANTILEVERED 
CONDITION 


Figures  6(i  and  67  present  the  V-g  plots  of  the  initial  analysis  of  the 
empennage.  Results  indicate  a  stable  system  throughout  the  proposed 
flight  envelope  for  anticipated  yaw  frequencies  of  10  and  20  cps.  Early 
preliminary  analyses  of  the  vertical  tail  with  and  without  the  effects  of 
an  elastic  horizontal  tail  showed  a  low  flutter  speed  which  necessitated 
increasing  the  torsional  stiffness  of  the  vertical  stabilizer.  Figure  9 
shows  the  increased  tip  stiffness  resulting  from  this  earlier  study.  The 
results  shown  in  Figures  66  and  67  reflect  this  latter  stiffness 
distribution. 

5.4  ANTISYMMETRIC  ANALYSIS  -  FREE- FREE  CONDITION 

Effects  of  fuselage  elasticity  and  aircraft  rigid  body  modes  on  the 
antisymmetric  flutter  characteristics  of  the  empennage,  as  shown  in 
Figures  68  and  69,  indicate  also  a  stable  system  throughout  the  design 
flight  envelope  of  the  XV-5A  aircraft. 

5.5  OVERALL  EVALUATION 

Basie  trends  shown  by  the  preliminary  empennage  flutter  analysis 
indicate  possible  existence  of  a  flutter  situation  which  is  dependent 
upon  establishment  of  the  actual  pitch  frequency  of  the  horizontal  tail. 
The  required  pitch  frequency  is  approximately  55  cps,  based  on  the 
results  of  the  symmetric  cantilevered  analysis.  Compressible  analy¬ 
ses  indicate  a  required  pitch  frequency  of  somewhat  greater  than 
55  cps;  however,  if  one  were  to  determine  the  required  pitch  frequency 
at  a  structural  damping  level  of  0.03  (see  Figures  59  through  61),  the 
required  frequency  would  be  of  the  order  of  55  cps.  In  all,  the  analy¬ 
ses  presented  in  this  report  should  be  taken  as  a  preliminary  evalua¬ 
tion  of  the  empennage.  Subsequent  analytical  and  experimental  investi¬ 
gations  should  be  directed  toward  finalization  of  the  actual  flutter 
characteristics  of  the  XV-5A  aircraft  with  due  regard  toward  the 
findings  presented  in  this  report. 
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Horizontal  Stabiliser  and  Elevator  ( Balanced)  -  Distribution  9A 


Node 

Weight 

Iba. 

Hi 

**  g 
Lb. -In. 

*Wv 

Ib.-In.2 

Ib.-In. 

«v 

lb.— la. 

0 

12.97 

972.6 

9541.9 

546.9 

17.61 

102.65 

1 

8.64 

81.0 

1848.0 

306.4 

10.08 

58.72 

2 

8.00 

70.0 

1966.1 

228.8 

6.52 

36.24 

S 

5.26 

67.0 

1070.6 

170.7 

5.20 

30.48 

4 

5.40 

51.0 

1008.5 

170.0 

6.31 

96.06 

6 

5.07 

42.6 

866.8 

145.8 

6.40 

97.00 

6 

9.00 

E31 

560.1 

03.6 

4.16 

24.94 

7 

9.09 

BStl 

584.8 

80.0 

3.71 

21.72 

8 

9.96 

It  "  o  ■■ 

958.7 

58.4 

2.06 

17.96 

8 

4.85 

470.9 

68.6 

/9.78 

22.16 

10 

■afl 

0.0 

0.0 

■D 

2 

♦58.77 

806.4 

11589.1 

1877.6 

66.67 

200.60 

Elevator  -  Distribution  9A 


♦Reflects  aa  elevator  wight  of  18,26  Iba.  including  balance  weights 
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MASS  AND  INERTIA  DISTRIBUTION,  HORIZONTAL  TAIL, 
ANTISYMMETRIC  ANALYSIS 


Horizontal  Stabilizer  and  Elevator  (Balanced)  -  Distribution  2 


Node 

Weight 

Lbs. 

Lb.  -In. 

Iv 

Lb.-In.^ 

•z  , 

IJb.  -In.2 

P 

uv 

Lb. -In.2 

su 

Lb. -In. 

Sv 

Lb. -In. 

0 

12.35 

299.0 

2907.8 

3006.4 

337.0 

14.52 

85.04 

1 

G.  97 

109.4 

2274.8 

2320.7 

243.7 

12.77 

74.81 

2 

6. 16 

76.9 

1669.8 

1700.6 

218.1 

8.40 

49.21 

3 

5.42 

61.1 

1160.1 

1184.6 

163.7 

5.73 

33.54 

4 

5.47 

45.1 

963.1 

983.9 

121.8 

5.68 

33.29 

5 

5.23 

39.3 

991.3 

1014.4 

104.3 

7.79 

45.64 

6 

4.07 

30.0 

618.9  ; 

635.5 

70.7 

4.83 

28.32 

7 

4.12 

25.7 

491.7 

505.4 

55.6 

4.02 

23.59 

8 

3.54 

25.1 

446.3 

459.9 

49.9 

3.96 

23.20 

9 

4.62 

44.9 

444.8 

469.4 

39.2 

3.96 

23.21 

10 

0 

2.2 

22.0 

23.2 

1.9 

0.20 

1.15 

2 

♦57.95 

758.7 

11990.6 

12304.0 

1405.9 

71.86 

421.00 

Elevator  -  Distribution  2 


Item 

Elevator  and  Balance  Weights 

Balance  Weights  Only 

Node 

Weight 

Lbs. 

2 

Lb.  -In. 

Sw 

Lb.  -In. 

Weight 

Lbs. 

Lb.  -In. 

Sw 

Lb. -In. 

0 

2.42 

30.79 

-0.37 

.544 

6.248 

-1.843 

1 

3.33 

62.50 

-1.96 

1.545 

17.552 

-5.204 

2 

2.14 

52.03 

-0.57 

.957 

9.593 

-3.029 

3 

1.98 

43.84 

-0.55 

.947 

8.602 

-2.852 

4 

2.10 

36.81 

-0.35 

.802 

6.495 

-2.282 

5 

2.25 

35.66 

-0.59 

.847 

6.302 

-2.310 

6 

1.81 

28.47 

-0.57 

.807 

5.261 

-2.059 

7 

1.87 

26.46  . 

-0.59 

.807 

4.595 

-1.924 

8 

1.85 

21.83 

-0.43 

.771 

3.706 

-1.689 

9 

.62 

5.40 

0.06 

.135 

0.591 

-0.283 

10 

0 

0 

0 

0 

0 

0 

2 

20.37 

343.79 

-5.92 

8.162 

68.945 

-23.475 

♦Reflects  an  elevator  weight  of  20.37  lbs.  including  balance  weights 
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MASS  AND  INERTIA  DISTRIBUTION, 
VERTICAL  TAIL,  ANTISYMMETRIC  ANALYSIS 


Vertical  Stabilizer  and  Rudder  (Balanced)  -  Distribution  1 


Node 

Weight 

Lb. 

^  2 
Lb.  -In. 

Iy  2 
Lb.  -In. 

p 

uv 

Lb. -In. 2 

su 

Lb.  -In. 

Sv 

Lb.  -In. 

0 

0 

0 

0 

0 

0 

0 

1 

0 

0 

0 

0 

0 

0 

2 

21.76 

981.0 

14796.1 

2661.7 

12.4 

77.39 

3 

16.89 

622.9 

9013.6 

1435.8 

2.1 

13.59 

4 

14.02 

437.1 

6164.9 

962.0 

4.1 

25.49 

5 

14.81 

405.3 

5394.1 

827.6 

2.3 

14.29 

6 

10.47 

265.3 

3283.8 

489.0 

0.3 

2.20 

7 

12.17 

264.0 

3442.0 

522.9 

0.3 

1.93 

8 

14.78 

297.2 

3913.3 

614.3 

-7.7 

-48.39 

9 

11.54 

281.0 

2783.6 

454.5 

-5.0 

-31.70 

10 

0 

0 

0 

0 

0 

0 

V 

w 

116.44* 

3553.8 

48791.4 

7967.8 

8.8 

54.80 

10 

102.3  t 

137926. 0t 

154726.0 t 

9621.0* 

92.53 t 

206. 26  t 

Rudder  -  Distribution  1 


Item 

Rudder  and  Balance  Weights 

Balance  Weights  Only 

Node 

Weight 

Lbs . 

*w  2 

Lb. -In. 

Lb.  -In.  ' 

Weights 

Lbs. 

*w  2 

Lb.  -In. 

sw 

Lb.  -In. 

0 

0 

0 

0 

0 

0 

0 

1 

0  (.0045) 

0.3 

-.0015 

0 

0 

0 

2 

6.66 

258.9 

.4284  ; 

i 

1.86 

54.47 

-10.34 

3 

3.78 

184.2 

-.  8654  , 

i 

2.39 

67.89 

-12.80 

4 

:  3.49 

171.3 

-.  8383 

1.83 

46.18 

-9.29 

5 

3.32 

111.4 

-.9953 

1.37 

30.10 

-6.41 

6 

2.06 

69.2 

-.6184 

.95 

18.20 

-4.16 

7 

2.46 

66.7 

-.6704 

1.13 

18.21 

-4.54 

8 

1.45 

33.9 

-.3380 

.59 

8.46 

-2.23 

9 

0 

0 

0 

0 

0 

0 

10 

0 

0 

0 

0 

0 

0 

2 

23.22 

895.9 

-3.8989 

10.12 

246.51 

-49.77 

♦Reflects  a  rudder  weight  of  20.  85  lbs.  including  balance  weights 

t  Reflects  horizontal  tail  lumped  mass  used  for  vertical  stabilizer 
modal  analysis 
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♦Treated  as  off -station  loads  in  fuselage  modal  analysis 


TABLE  5  (Continued) 


MASS  AND  INERTIA  DISTRIBUTION,  FUSELAGE, 
ANTISYMMETRIC  ANALYSIS 


Item 

Total  Fuselage 

Node 

Weight 

Lbs. 

Ix  2 

Lb. -In. 2 

h  2 

Lb. -in. 

SX 

Lb.  -In. 

0 

14.30 

484 

1 

80.25 

34018 

2 

188.75 

84267 

3 

287.47 

126017 

4 

877.18 

387156 

5 

888.41 

269249 

6 

1025.23 

434862 

7 

733.38 

438267 

8 

709.95 

345643 

9 

860.91 

251607 

10 

287.40 

128358 

11 

121.80 

58261 

12 

93.48 

31011 

13 

44.97 

11955 

14 

26.36 

3442 

15 

7.67 

253 

1 

* 

1238.27 

6479172 

7184964 

-14188 

* 

1438.00 

7528472 

8348565  j 

-16486 

* 

218.74 

1342264 

445814 

i 

2  * 

9142.52 

17954758 

15979343 

I 

-30674 

♦Treated  as  off-station  loads  in  fuselage  modal  analysis 


TABLE  6 


MODAL  DATA,  HORIZONTAL  STABILIZER,  SYMMETRIC  ANALYSIS 


B.L. 

Node 

1st  Mode 
f  =  29.5  cps 

2nd  Mode 
f  =  91.4  cps 

3rd  Mode 
f  =  110.9  cps 

<P 

6x  102 

<P 

© 

<P 

© 

3.03 

0 

0 

0 

0 

o 

0 

0 

10.637 

1 

.009 

.0336 

-.0005 

-.0072 

-.0524 

-.00445 

18.244 

2 

.036 

.0733 

.0015 

-.0154 

-.1965 

-.00892 

25.851 

3 

.083 

.1204 

.0108 

-.0245 

-.4013 

-.01264 

33.458 

4 

.152 

.1773 

.0379 

-.0345 

-.6223 

-.01448 

41.065 

5 

.244 

.2441 

.0953 

-.0449 

-.7930 

-.01179 

48.672 

6 

.360 

.3185 

.1952 

-. 0552 

-.  8305 

-.00131 

56.279 

7 

.500 

.4249 

.3463 

-.0685 

-.6581 

.02319 

63.886 

8 

.658 

.5486 

.5447 

-.0824 

-.2460 

. 06244 

71.493 

9 

.828 

.6533 

.7701 

-.0936 

.3506 

.  10275 

79.100 

10 

1.000 

.  6533 

1.0000 

-.0936 

1.0000 

. 10275 
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TABLE  7 


MODAL  DATA,  FUSELAGE, 
SYMMETRIC  ANALYSIS 


F.  S. 

Node 

1st  Mode 
f  =  8. 2  cps 

2nd  Mode 
f  =  15. 2  cps 

3rd  Mode 
f  =  23. 2  cps 

4th  Mode 
f  =  33.0  cps 

<P 

<P 

<P 

<P 

0 

0 

1.0000 

1.0000 

1.0000 

1.0000 

34.33 

1 

.7260 

.5474 

.3354 

.0745 

68.66 

2 

.4955 

.2433 

.0202 

-.1329 

102.99 

3 

.2890 

.0120 

-.1548 

-.1440 

137.33 

4 

.1211 

-.1146 

-.1637 

-.0191 

171.66 

5 

.0035 

-.1222 

-.0420 

.0764 

205.99 

6 

-.0794 

-.0767 

.  0735 

.0745 

240.33 

7 

-.1164 

-.0053 

.0987 

-.0028 

258.34 

< 

8 

-.1167 

.0288 

.0725 

-.0336 

274.66 

9 

-.  1043 

.0535 

.0247 

-.0409 

308.99 

10 

-.0498 

.0860 

-.0851 

-.0186 

343.33 

11 

.0289 

.0992 

-.1943 

.0389 

377.66 

12 

.1262 

.0913 

-.2695 

.1000 

411.99 

13 

.2543 

.0464 

-.2690 

.1414 

446.33 

14 

.4083 

-.0329 

-.1829 

.1426 

459.37 

15 

.4739 

-.0735 

-.1204 

.1249 

480.66 

16 

.5789 

-.1396 

-.0062 

.0780 

514.99 

17 

.7539 

-.2535 

.1927 

.0002 

TABLE  8 


MODAL  DATA,  HORIZONTAL  STABILIZER, 
ANTISYMMETRIC  ANALYSIS 


B.L. 

Node 

1st  Mode 
f  =  29.4  cps 

2nd  Mode 
f  =  90 . 2  cps 

3rd  Mode 
f  =  111.0  cps 

<P 

ex  io2 

<P 

e 

<P 

© 

3.03 

0 

0 

0 

0 

0 

0 

0 

10.637 

1 

.0089 

.0368 

-.0021 

-.0068 

-.0694 

-. 0045 

18.244 

2 

.0365 

.0805 

-.0049 

-.  0144 

-.2613 

-.0088 

25.851 

3 

.0838 

.1322 

-.0016 

-.0226 

-.5361 

-.0117 

33.458 

4 

.1526 

.1947 

.0180 

-.0314 

-.  8357 

-.0118 

41.065 

5 

.2450 

.2692 

.0679 

-.0407 

-1.0747 

-. 0063 

48.672 

6 

.3616 

.3526 

.1625 

-.0493 

-1.1495 

.0100 

56.279 

7 

.5012 

.4721 

.3121 

-.0597 

-.9632 

.0463 

63.886 

8 

.6594 

.6141 

.5153 

-.0700 

-.4771 

.1041 

71.493 

9 

.8285 

.7255 

.7536 

-.0768 

.2318 

.1585 

79.100 

10 

1.0000 

.7320 

1.0000 

-.0771 

1.0000 

.1617 
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FUSELAGE,  ANTISYMMETRIC  ANALYSIS 
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Figure  6  Planform  Horizontal  Tail 


Figure  7  Planform  Vertical  Tail 
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Figure  9  Stiffness  Distribution  -  Vertical  Stabilizer 
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Figure  10  Stiffness  Distribution  -  Fuselage 
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Figure  13  Mode  Shapes  -  Horizontal  Stabilizer  -  Symmetric  Analysis 
(Mode  3  -  £  ■  110. 9  cps) 


Figure  14  Mode  Shapes  -  Fuselage  -  Symmetric  Analysis 
(Mode  1|  f  "  8. 2  cps  •  Mode  2,  t  •  15*  2  cps) 
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Figure  21  Mode  Shapes  -  Vertical  Stabilizer  -  Antisymmetric  Analysis 
(Mode  3  -  f  =  27. 9  cps) 
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Figure  22  Mode  Shapes  -  Fuselage  -  Antisymmetric  Analysis 
(Mode  1,  f  =  5.6  cps  -  Mode  2,  f  =  13.6  cps) 
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Figure  23  Mode  Shapes  -  Fuselage  -  Antisymmetric  Analysis 
(Mode  3,  f  =  20.6  cps  -  Mode  4,  f  =  26.  2  cps) 
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Figure  24  Mode  Shapes  -  Fuselage  -  Antisymmetric  Analysis 
(Mode  1,  f  =  13.9  cps  -  Mode  2,  f  =  19.0  cps) 
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Figure  26  Symmetric  Analysis  -  Cantilevered  Condition 
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Figure  28  Symmetric  Analysis  -  Cantilevered  Condition 


Figure  31  Symmetric  Analysis  -  Cantilevered  Condition 
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Figure  32  Symmetric  Analysis  -  Cantilevered  Condition 
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Figure  34  Symmetric  Analysis  -  Cantilevered  Condition 
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Figure  36  Symmetric  Analysis  -  Cantilevered  Condition 
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Figure  37  Symmetric  Analysis  -  Cantilevered  Condition 
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Figure  39  Symmetric  Analysis  -  Cantilevered  Condition 
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Figure  44  Symmetric  Analysis  -  Cantilevered  Condition  -  Parameter  Effects 
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Figure  45  Symmetric  Analysis  -  Cantilevered  Condition  -  Modal  Coupling 
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Figure  49  Symmetric  Analysis  -  Cantilevered  Condition 
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Figure  53  Symmetric  Analysis  -  Cantilevered  Condition 
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Figure  60  Symmetric  Analysis  -  Free-Free  Condition 
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Figure  62  Symmetric  Analysis  -  Free-Free  Condition 
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Figure  63  Symmetric  Analysis  -  Free-Free  Condition 
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Figure  64  Symmetric  Analysis  - 
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Figure  66  Antisymmetric  Analysis  -  Cantilevered  Condition 
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Figure  67  Antisymmetric  Analysis  -  Cantilevered  Condition 
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INCOMPRESSIBLE  _ 
SEA-LEVEL 
YAW  -  10  CPS 
RUD,  -  to  CPS 
ELEV.  -  15  CPS 
FLAP  -  10  CPS 


VELOCITY  -  KNOTS 


Figure  68  Antisymmetric  Analysis  -  Free-Free  Condition 
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